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A PROGRAM  FOR  STORING  OCEANOGRAPHIC  DATA 
ON  MAGNETIC  TAPE 


1.0  IDENTIFICATION 

1.1  Title 

Storage  Program  for  Navigation,  Bathymetric,  and  Magnetics  Data  on 
Magnetic  Tape. 

1.2  Identification  Name 
GEODATA. 

1.3  Classification  Code 
None. 

14  NRL  Research  Computation  Center  Identification  Nnmhor 
None. 

1-5  Entry  Points 
GEODATA. 

1-6  Programming  Language 

Language:  CDC  3600/3800  Fortran. 

Routine  Type:  Program. 

Operating  System:  Drum  Scope  2.1. 

1-7  Computer  and  Configuration 

CDC  3800. 

1.8  Contributor  or  Programmer 

Marilyn  L.  Blodgett,  Code  4223MB,  Researc  * Computation  Center, 
written  for  Environmental  Sciences  Section,  Acoustics  Division. 

1.9  NRL  — Naval  Research  Laboratory,  Washington,  D.C.  .20375. 


NOTE-  Manuscript  submitted  January  9 1975. 
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1.10  Program  Availability 

If  supplied  with  a magnetic  tape,  a copy  of  this  program  will  be  made 
available  by  the  Environmental  Sciences  Section,  Acoustics  Division. 

1.11  Verification 

This  program  has  been  used  and  tested  by  the  Environmental  Sciences 
Section,  Acoustics  Division,  for  the  past  year. 

1.12  Date 
July  1974. 

2.0  PURPOSE 

2.1  Description  of  the  Routine 

This  program  transfers  the  data  collected  by  an  oceanographic  or  geo- 
physical experiment  to  magnetic  tape.  The  program  will  take  up  to 
3000  navigational  data  points  and  an  unlimited  number  of  bathymetric 
and  magnetic  data  points.  The  program  is  set  up  so  that  if  one  tape  is 
filled,  any  number  of  continuation  tapes  can  be  used.  The  data  tape 
v.iU  have  one  logical  record  (of  80  characters)  for  each  data  point.  The 
different  types  of  data  (navigational,  bathymetric,  and  magnetics)  will 
be  separated  by  an  end-of-file  mark,  with  a double  end-of-file  mark  at 
the  end  of  all  the  data. 

2.1.1  Navigation  Data 

This  program  is  presently  set  up  to  read  the  navigation  data 
from  cards.  The  data  are  in  degrees,  minutes,  and  hundredths 
of  a minute;  each  navigation  fix  has  an  associated  Julian  date 
and  time  (24-hour  clock).  The  southern  latitudes  and  the 
western  longitudes  are  preceded  by  a minus  sign.  Thi  pro- 
gram changes  the  Julian  date  to  month  and  day  and  changes 
the  minutes  and  hundredths  of  minutes  of  latitude  and  long- 
itude to  ter  thousandths  of  a degree. 

2.1.2  Bathymetry  Data 


The  bathymetry  data  are  read  in  with  the  Julian  date  'ollowed 
by  five  time-depth  groups  to  a card.  Time  is  read  to  tenths 
of  a minute.  The  depth  can  be  either  uncorrected  fathoms 
or  uncorrected  meters,  but  not  both  in  the  same  run.  This 
program  converts  the  Julian  date  to  month  and  prints  out 
uncorrected  fathoms,  unconnected  meters,  corrected  meters, 
and  Matthews  zone.  The  uncorrected  meters  arc  not  written 
on  the  data  tape. 
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2.1.3  Magnetics  Data 

The  magnetics  data  are  read  in  with  the  Julian  date,  then  the 
hour,  followed  by  12  magnetic  data  points  (one  every  five 
minutes).  The  program  converts  the  Julian  date  to  month  and 
day  and  converts  the  total  magnetic  intensity  in  gammas  (the 
data  read  in)  to  residual  magnetic  intensity  using  the  Inter- 
national Geomagnetic  Reference  Field  formula. 

2.1.4  Gravity  Data 

This  program  is  not  set  u*'  to  take  gravity  data,  but  a gravity 
subroutine  could  easily  be  'rted. 

2.2  Problem  Background 

The  rapid  accumulation  of  many  boxes  of  oceanographic  data  nn  com- 
puter cards  necessitated  the  transfer  of  the  data  from  cards  to  magnetic 
tape.  It  was  decided  that  the  format  recommended  by  the  National 
Research  Council  of  the  National  Academy  of  Sciences  would  be  used 
with  or."  slight  modification,  the  addition  of  a fix  number  with  the 
navigation. 


3.3.1  Unique  Storage 

36533  Octal  (15707  decimal)  locations  exclusive  of  system 
library  functions. 

3.3.2  Common  Blocks 
None. 

3.3.3  Temporary  Storage 
None. 

3.4  Messages  and  Instructions  to  the  Operator 
None. 
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3 5 Error  Returns,  Messages,  and  Codes 

Deck  set  up  incorrectly. 

More  than  3000  navigation  data  cards. 

Unable  to  interj>olate  latitude  and  longitude. 

3-6  Informative  Message:  to  the  User 

None. 

3.7  Input 

Input  parameters  and  data  are  read  in  via  punched  cards.  The  format 
statements  can  easily  be  changed  to  meet  requirements  of  the  user, 
bee  Appendix  A for  samples  of  our  data  formats.  Each  data  type 
(navigation,  bathymetry,  magnetics)  are  bracketed  by  a control  card 
with  a negative  number  in  columns  1 and  2.  The  control  card  preceding 
the  uata  sets  up  the  different  options  available.  Appendix  B is  a com- 
plete description  of  the  input  setup. 


3.8  Output 

The  program  prints  out  on  the  standard  printer  (logical  unit  61)  the 
ship  and  cruise  identification,  number  of  cards  read,  number  of  logical 
records  written  for  navigation,  bathymetry,  and  magnetic  data.  There 
is  an  option  for  listing  all  the  records  written  on  the  output  tape  and 
for  punching  out  cards.  Appendix  C shows  samples  of  the  output. 

3.9  Formats 

Appendix  B shows  the  program  deck  structure. 

3-10  External  Routines  and  Symbols 

SKIPt-TLE,  BACXFILE,  XM0DF,  SQRTF,  SINF,  CQSF. 

3.11  Timing 

This  program  will  process  approximately  1300  data  points  a minu,e  If 
a listing  of  all  processed  data  is  required,  each  data  point  will  generate 

Conusec*u?nt,y-  **  Program  can  process  a great  deal  of  data 
in  little  time  but  also  produce  large  printouts. 

3.12  Accuracy 

Not  applicable. 

3.13  Cautions  to  Users 
None. 
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3.14  Program  Deck  Structure 
See  Appendix  B. 

3.15  References  — Literature 

“Formats  for  Marine  Geophysical  Data  Exchange,”  National  Academy 
of  Sciences,  June,  1972. 

METHOD  OR  ALGORITHM 

Not  applicable. 

FLOW  CHART  AND/OR  SOURCE  LANGUAGE  LISTING 
Flow  chart  and  listing  are  given  in  Appendix  D. 

COMPARISON 

There  are  no  other  known  programs  available  for  comparison. 

TEST  METHOD  AND  RESULTS 

A sample  of  the  three  types  of  records  written  on  the  output  tape  are  included 
in  Appendix  C. 

REMARKS 

The  authors  thank  Gary  Flenner  for  his  review  of  this  report  and  Wayne  Worsley 
for  preparing  the  illustrations  in  the  appendixes.  The  authors  also  thank  Leon 
LaLumiere  for  providing  some  of  the  subroutines  used  by  this  program. 
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Magnetics  Input  Format 


302  READ  (60,104)  IFN  , JUDY , RMIIR , ( IMG  ( I ) , 1 = 1 , 1 2) 
304  FORMAT  (12, 13, 1XF2, 12(1X15)) 


Appendix  B 

DECK  ASSEMBLY 


16  End  of  file  card 
1 5 Finish  card 

Magnetics  end  card 


Magnetics  data  deck 
Magnetics  parameter  card 


^ ^ 9B 


Bathymetry  end  card 


Bathymetry  data  deck 


“ 7S 

Y J 

* 9A 

- 

10 


Zone  change  card 


' 9A  Bathymetry  data  deck  (portion  of 

data  taken  at  new  Matthews  Zone) 

J-2.  Zone  change  card  (needed  only  if  change 
in  Matthews  Zone) 

Bathymetry  data  deck  (portion  of  data  deck 
taken  at  Matthews  Zone  specified  on  bathymetry 
Bathymetry  parameter  card  parameter  card) 


I 7 Navigation  end  card 

T 

S 6 Navigation  data  deck 

5 Navigation  parameter  card 


Run  card 


Program  geodata  (binary  deck)* 


Equip  card 


Job  card 


*If  the  FORTRAN  source  eck  is  used  instead  of  the  binary  deck,  a 
FORTRAN  card  is  require  after  the  Equip  card.  The  FORTRAN  card 
reads--7/9  FTN,  L,  R,  X.  In  addition  a SCOPE  card  with  SCOPE  starting 
in  column  10  and  a LOAD  card  which  reads  7/9  LOAD  must  follow  the 
source  deck. 
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Number 


Card  Column 

Title  Number  Description 


1 


2 


3 


4 


5 


Job  1-21 


Equip  1-18 


Program  Deck  of 

GEODATA  cards 


Run  1-13 


Navigation  1-2 

Parameter 


5 


7/9  JOB,  charge  No.,  ID  No.,  time. 
The  charge  number  consists  of 
alphanumeric  characters  of  any 
length.  The  ID  number  is  the  pro- 
grammer identification;  it  may  be 
any  length  and  appears  as  given  in 
the  control  card  listing.  The  time 
is  the  maximum  time  limit  in  min- 
utes for  the  entire  job.  See  page 
2-2  of  the  3600/3800  Computer 
Systems  Drum  Scope  Manual. 

7/9  EQUIP,  10  - **,  WO,  HI 
10  = logical  unit  number. 

WO  = write  only. 

HI  = high  density. 

This  is  the  main  program  with  asso- 
ciated subroutines.  If  the  Fortran 
source  deck  is  used  instead  of  the 
binary  deck,  a Fortran  card  is  re- 
quiret  after  the  Equip  card.  The 
Fortran  card  reads  7/9  FTN,  L,  R, 

X.  In  addition  a Scope  card  with 
SCOPE  starting  in  column  10  and  a 
LOAD  must  follow  the  source  deck. 

7/9  RUN,  T,  P,  R,  M,  D 
T = time  lirmt  in  minutes. 

P = maximum  number  of  print  or 
write  operations.  R,  M,  D may  be 
left  blank.  See  page  2-15  of  the 
3600/3800  Computer  Systems  Drum 
Scope  Manual. 

-1 

This  number  designates  the  Naviga- 
tion Parameter  card. 

0,  1,  2 

0 = write  the  navigation  data  on  a 
new  tape. 

1 = read  the  navigatirn  data  cards 
and  store  the  information  necessary 
for  interpolation.  The  program  will 
not  write  the  navigation  data  on  the 
tape.  It  will  only  write  the  bathy- 
metry data  on  a tape  that  already 
contains  the  navigation  and  magnetic 
data. 
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Card  Column 

Number  Title  Number 


19-20 


23-30 


6 Navigation 
Data  Deck 

7 Navigation  1-2 

End 

8 Bathymetr'  1-2 

Parameter 

4-5 
14  io 

11 


Description 

2 = Read  the  navigation  data  cards 
and  store  the  information  necessary 
for  interpolation.  The  program  will 
not  write  the  navigation  data  on  the 
tape.  It  will  only  write  the  mag- 
netics data  on  a tape  that  already 
contains  the  navigation  and  bathy- 
metry data. 

0,  1,  2,  3 

0 = just  write  the  navigation  data 
on  tape  in  the  required  format. 

1 = write  the  navigation  data  on 
tape  and  print  out  a listing  of  the 
data. 

2 = write  the  navigation  data  on 
tape  and  punch  out  a card  for  each 
logical  record  written  on  the  tape. 

3 = write  the  navigation  data  on 
tape  and  print  out  a listing.  In  ad- 
dition, punch  out  a card  for  each 
logical  record  written  on  the  tape. 

If  card  columns  4-5  are  1 or  2,  more 
of  the  options  in  card  columns  19-20 
are  available. 

Columns  23-30  are  for  cruise  identi- 
fication. Any  eight-digit  alphanu- 
meric number  may  be  used. 

Place  the  Navigation  Data  Deck  after 
the  Navigation  Parameter  card. 

-5 

This  number  designates  the  end  of 
the  Navigation  Data  Deck. 

-2 

This  number  designates  the  Bathy- 
metry Parameter  card. 

Columns  4-5  are  for  the  Matthews 
zone. 

1 or  2 

1 = Depth  data  in  uncorrected  fathoms. 
r.  = Depth  data  in  uncorrected  meters. 
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Number 


9 


10 


11 


12 


Card  Column 

Title  Number 


19-20 


Bathymetry 
Data  Deck 


Zone  Change  1-2 


4-6 

Bathymetry  1-2 

End 

Magnetics  1-2 

Parameter 

6-10 


Description 


0,  1,  2,  or  3 

0 * just  write  the  bathymetry  data 
on  tape  in  the  required  format. 

1 * write  the  bathymetry  data  on 
tap''  and  print  ou<  a listing  of  the 
data. 

2 = write  the  bathymetry  data  on 
tape  and  punch  out  a card  for  each 
logical  record  written  on  the  tape. 

3 = write  the  bathymetry  data  on 
tape  and  print  out  a listing.  In  ad- 
dition punch  out  a card  for  each 
logical  record  written  on  the  tape. 

Place  that  portion  of  the  Bathy- 
metry Data  Deck  which  was  taken 
at  the  Matthews  zone  specified  on 
the  Bathymetry  Parameter  card  after 
the  Bathymetry  Parameter  card. 

—6 

This  number  designates  a change  in 
the  Matthews  zone.  A Zone  Change 
card  must  precede  the  first  data  card 
taken  at  the  new  Matthews  zone. 
These  Zone  Change  cards  can  be 
scattered  throughout  the  Bathymetry 
Data  Deck. 

Value  of  the  new  Matthews  zone. 

-7 

This  number  designates  the  end  of 
the  Bathymetry  Data  Deck. 

-3 

This  number  designates  the  Mag- 
netics Parameter  card. 

Height  in  feet  above  or  below  mean 
sea  level. 


19-20  0,  1,  2,  or  3 

0 = just  write  the  magnetics  data 
on  tape  in  the  required  format. 

1 = write  the  magnetics  data  on 
tape  and  print  out  a listing  of  the 
data. 
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Number 

Card 

Tide 

Column 

Number 

Description 

2 = write  the  magnetics  data  on 
tape  and  punch  out  a card  for  each 
logical  record. 

3 = write  the  magnetics  data  on 
tape  and  print  out  a listing.  In  ad- 
dition, punch  out  a card  for  each 
logical  record  written  on  the  tape. 

13 

Magnetics 
Data  Deck 

Place  the  magnetics  data  deck  after 
the  magnetics  parameter  card. 

14 

Magnetics 

End 

1-2 

-8 

This  number  designates  the  end  of 
the  Magnetics  Data  Deck. 

15 

Finish 

1-2 

— 4 

This  number  designates  the  end  of 
all  data. 

16 

End-of-File 

■ ' • -- 
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M X OOOOOOOOO  O w4  w4  vi  H H M (V  <V  N CV  ^ ^ J 

u.z3t^<vrx<vrM^rwr«tvf«wrwrgrw<xfvrgrgfv«>#cxfx<x»xcvr.  o#cx«xcxc.rgrw<xcxcxcx«N 


t m O C»  C OOO  o o #o  O » c » 


:o  c cff  c cc  c c c co< 


:cccccccc < 
k KMIMr  «\  tfMIMIMn 


O Opv  «MO  (V« 


— MrtMNN  HO 


n ok  n n om  oo  o 


)««••••••  ,'•••••••••• 


_j«rvr*r*r*r*r*r*r*««'«c:  c 


fs.  oonn^cn^i 
•>  hit  •>  n * m « • « 


4 NWNNNNNNM 


»r»  fw  ir  « »o  c.  or 

OMTK  N • • ^ p>  OlVv'** 

r>.ON'0pvP.|r>0««M»«  N »ON  ^ 

w«wnr)w«Kiw«i 


««• o«oc 

I « • O HO  O 
I H H MK]  « O O 


M «oooo«oo«r»*«  °5»  00*^000 

x n«uv«  hio  «<n  n » * h«  <v  »nr  *««h«  n <*  **  •*  * 


OOOHHHNMI 


OOHHNN«W«OOH( 

IHHHHHHHHHHHHHHH  H N (V  N N N NN  N N 


NNNNWWNNNNNNNNNNNNNNNNNNNNNNNNNNWNPlrtMWOO 

NNM  l«N  rin  N N MM  N(4I(«  111  Min  M M W W U»  W W W •*  W M M W M M "•  1N*  w ,V  w ,V 
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u u r.  < a A r « r • » « < c c ■ r c c v v 
**-**'«*■>»  o r>.  iv  n « ir  « ^ ■ » o oo 

“ » (vfvrvNKrvrvrv.vrvrvfvNcvevcvMvfvfv 


r r < r n m * m • 

HHtvNfoarc^  o hmn  n t f » v a « « ^ v a kn  < 


uw 

**  • c c ■*  c c c c •«  c c 
out  • • • • -*•■■•« 


* «■  *cc  cc  c c c c.cc  c«c«c  e<c  c c c.c  c cm 


► v. 

U X 

icbc«*b<o«-»  ^ r-  r-  r\  b-  .-#  as 
B 

o ► o«  n • b*  p"  • .*  r-  b,  r-  *n 


H(MNMNnr>rjr)«  » 000 

•t'  r r'  r-  r-  rv  «■  ■•  rrrr  rrrr  rr  r « « « 


4CC^-fV<««  rc^rvrvMvca^mivt 
**  ^ « n n <»  r c n « * a a nn  o»  » « b-  < 

coooooooHHNnnnnnnNNNNt 

<««««««« 


sSsxsssssS«s»ss&::£*fi&«£cs;:3*srssss;s:8f»s?s;ss;;: 


i c cic’  ccctctocc  Cjor  cooccccoccoc.bococoB»B»( 

-*  *****•••••••••••••»••••••••••«••• 


3 •«  ••«  r c r a« 


ubub  r • « • n « a ^ ••  «••«  •«•  «»  « 

' *•  «*•  • A*  Br  n « « MB  « BB  « PIB  B I 

i^^bbbbbbbbbb nH-maap • b • •* i 


: ??  ;s  « •?  ?s?t  °r  ••isisss  w?s-«is 

:ssasasaa533533iara5TO5ratiti,ai««^5^wm«m«'6?i'a?P 


u (X(Mo.o,ni(v<o*(xra(^ra(xrao«(v. 


r>p>p>r>p>P'p.p.p. 


«f DonooooDoonDOoooonooononoooooooooooooo 


oaonnnnooonoDooonoooo 


*j-*r><rv,rvf><rxrNfNrvr'.rvr>«f\<*>r.rv.rv»vrv»>*,x*vf 

l.««r.o«0f,„„f,r,c,t,0tJ0CJt|(||1(1( 


**•"  0 cn?4 
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to1 

Z ► 


itnr*«rh«  « ntv  « o«  oc  iv  « r>  rw  o r»  «*  « « o o po  r-  a.*:-  4 r- » » n « » h ► Ho«ro<eti*  »•»  o < c * c «n<r  if  *•' 


o on««  «4«  »«s  *4«  « ok  »<*NHHHoo«HHHoo*»>.r'»nno^«  o r-  4 « nrMHucsrortv 

3 • r>4«orv»r'4  HO*h«vN*,iNO«««nf«ot«HCH<KO««,ioM>)c«iMK«  *<►«  *<r««\r.  ^*Nrirrirr*w 

*-  4r»nwH»«r>4tfMnNHO4«^r*nMHO4r'<i«'0HO«N«in>>ir»rto«>«r«NH3cr-«irK)ivHO©««^H 

— o-o  o-o  o oo  oooooooooooo®oo«r'r-r*r.r*»*»**-oo«o«*««ooooo 


► 4 o ® or»*»n*r^o 

w»  r *<•  » HN  (0  4 4 

Kr  4f^r>  ••  4 o ««  •«< 


ovi  H4nvnP-(v4  w*  <r  04  « cv  4 h»niNNM«M< 
»«m»«f400f*®o*4«*or*r>»o<M»«’oor*« 
««»44N4«400HrwNn4«  «\«N  4M>C»  OC-i 


h Hfyr>,(VNrn«  r « 


0000000000 


000000000004%  00  00  00  00  OO  0000000000  00  OOOO  oooooooooooooooo  00  000 
40®  h hmcmp)  n ••  « 00  »*»<nMnn««r>>  ^*«cviMnr)«  «r  r ^ hiv  jvrxo  « 400 


^u*  p^  p^  p^  p^  p"*  p^  p"*  r*  p1*  p1*  p^  p^  p^  p^  P1*  r»  p^  p^  p^  p*  p^  p**  p*  p**  p^  p^  p*.  4*  p^  p^  **  p^  p*  p*  p^  p^  p*  p^  p*  p*  *■•  P1*-  p*»  P1*-  p*»  p*  p**  p*»  p**  p**-  r*  p»  p«*  **■  p**  p* 
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Sample  Printout  of  Data  Processed 


SWIP  AN”1  CRUISE  IDENTIFICATION  7^1602 

NAVIGATION  DATA 

number  navigation  cards  read  p &e 

number  logical  records  written  vn  tape  = 60 

BaTHTMEtRV  data 

DEPTH  D.'Ta  GIVEN  IN  UNCGRRECTED  MfTFRS 
MATTHEW^  7ENES  PASSED  THROUGH  WFRE  123 

number  ^f  bathymetry  cards  read  p 7* 

NUMBER  -F  LOGICAL  RECORDS  WRITTEN  « 357 

MAGNETKS  CAT* 

HEIGHT  ; N FFET  ABO'.E  OR  BELOW  SFA  LEVFL  IS  *20 
NUMBER  F MAGNETICS  CARDS  READ  * ''3 

NUMBER  °F  LOGICAL  RECORDS  WRITTEN  « 258 


18 


Number 


229  WRITE  (62,224)  (CRUISE  , ITMZNE  , IYEAR  ( I ) ,MONT!I(I)  , I DAY  (I)  ,IIR(I)  ,XMIN(I 
1)  , XLAT  ( I ) , XLONG  ( I ) ,NAVFIX(I)  ,ICRMET(I) ,IFIX(I) ,I=1,IWRITE) 

224  FORMAT  (A8 , 15 , 1 2 , 1 2 , 1 2 , IXF2 , F3 ,F8 . 4 , F9 . 4 , 10X1 5 , 1 5 , 1 2 , 16X) 


4 

Implies  a decimal  point 


I 


20 


5 5 i 5 5 5 5 
F66R6|56 

Minn; 


99191991 

- M •>  » 
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Magnetics  Output  Card 


a o o o o o l 

i,i«  < 

1 1 1 i 1 1 1 

1112117 

mini 


a o a | q a a o 


32  3 WRITE (62, 318)  (CUR1SE , IT'tZNE , I YEAR ( I ) , MONTH (I ) ,IDAY(I)  ,HR(I)  , XMIN (I 
1)  , XLAT ( I ) , XLONG ( I ) ,NAVFIX(I)  , ICRMET ( I ) , I=1,IWRITE) 

318  FORMAT  (A8 ,1  5, 12 , 12 ,12 , 1XF2  ,F3  ,F8 . 4 ,F9 . 4 , 23X  , 15 , 15 , 5X) 


Implies  a decimal  point 


A 


Appendix  D1 
FLOW  CHART 


22 


i 


i 


j 


rt  r»  r»  o n n o n noon 


Appendix  D2 


SOURCE  LANGUAGE  LISTING 


PROGRAM  GEODATA 

PROGRAMMER  MARILYN  L.  BLODGETT  CODE  7817MLB 
NAVAL  RESEARCH  LABORATORY 
WASHINGTON  p.C.  20375 

DIMENSION  IYEAR(20)|MeNTH(20)»lDAY<20)‘,HR(20),XMIN(26)',XLAT(20'). 
lXL0NG(2o),iriX(2O),NAVFlx(2O).lDATE(3ooO)»STlME(3OOO).STL'At(3O6O). 
2STL0NG<3OOO>.KEEPMTI24)»BMR(5),BM1N{5)'.  I DP  (5 ) . ICRM6T 120  > , ZAVF IX  C 20 
3 > • I TH ( 12 ) | I MG  < 12 ) . ZH I N < 20  > , SI  DPI  20 5 
CATa( (I TM| I ), I ■1.121 «0.3.l0 1 15,20,25 1 30, 35,4 0.45,5 0.55) 

REAL  IDEPTH 
REWIND  06 
REWIND  10 
IFLIGMTil 
ITMZNE'O 
ICNTBiO 
1CNTN*0 
1CTmAQ«0 

READ  A PARAMETER  CARD 

160  READ<6O.lO)IC0DE.MT|AlTt INPUT,I0UT,SHJPID 
1C  FORMAT! 1 2, 1X12. F5, 3X12,3x1 2, 2XA0) 

I F < I CODE 112.11.11 

11  WR I TE ( 6l , 13 ) 

13  FORMAT  (1H0.23HDECK  SET  UP  INCORRECTLY) 

ST0P 

12  IF(4*IC0DE)11,14.14 

14  ICODEq - 1 CODE 

GO  TO  ( 100,200,300.400)  1 I CODE 

PROCESS  THE  NAVIGATION  DATA  CARDS 

ICO  CRUISE'SHIPID 
I0PT»IOuT 
NOWRlT'MT 

IF(N0WR1T  ) 166 ,165(166 
ie«  CALL  SKIPFILE(IO) 

READ(10,110)CRUISE 
IF (EOF  ,10  3167 , 166 
167  CALL  BACKFILE(IO) 

GO  TO  120 

165  IF< I OPT  >121.120 ,121 

121  I F < 10PT-2)122,120,122 

122  wR I TE ( 6i , 3 01 ) 

5 C 1 F6RMAT(lHl,i23HSHJP  AND  CRUISE  TIME  YEAR  MONtH  DAY 

1 HOuR  MINUTE  latitude  longitude  f I X 

2 Fix) 

WRITE (61,502) 

5c2  FORMAT ( lH  ,24HIDENTIFICAtION  Z0NE , 79X , 22HDESCR ( PT I ON  NUMB 

1ER5 

120  IENDFL'O 
I WR  I TE*0 
NLR£C«0 
NPREC^O 
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r 


NRL  REPORT  7861 


1ST0R6-1 

C 

0 READ  A NAVIGATION  CARD 

ice  read(60,iodisn,iyi vudyj rnmr  , rnmin, i r n*» xLap* xlam, xlod, xlom  , i r i x r y 

ClCl  FORMAT  (12, 4X 12. 1XI3|1XF2,F2, 2X13.  lXF2;F6.2,F4,r6,2‘,  1X12. 36X) 

1C1  FERMAT  (12,4X12.  IX 1 3I1XF2.F2. 15, F3.F6,2'.F«.F6.2‘.  1X1 2) 
IF(5«1SN)11,102,104 
1 C 2 IF(NCHR1T)160,170,160 

170  1ENDFU«1 

1F(1WR1tE)150, 150,103 

103  NlREC»NtREC*tWRlTE 
NPREC«NPREC^1 
1BEG*Iwr1TE*1 

D0  117  k«IBEG,20 
1 rE AR (K ) «0 
MeNTH(K)»0 
IDAy  CK » *0 
HR  ( K ) fl 0 
XMlN(K)«0 ,0 
2M I N (K » *0 
XLAT  CK » *0 . 0 
XL0nG(K)«O.O 
IF  1 X (K ) *0 
117  N AVF I X ( K ) <0 
1 WR 1 TE«20 
G0  TO  140 

104  1 YR» I Y 
1CNTN^ICNTN*1 
1F(3001-ICNTN)105,105U06 

105  WR!TE(61,107) 

107  F0RMAT(lMi,36HM0RE  THAN  3000  NAVIGATION  DATA  CARDS) 

STOP 

C 

C COMPUTE  DAY  and  MONTH  FROM  JULIAN  DATE 
r 

ice  CALL  JUlIAN(  IYR,  JUDY. lOl IM.LRYR) 

C XL0D««XlOD 
XLOD««XlOD 

MlNuS*4oOOOOOOOOOOODOOB 

XLAm«A0SF(XLAM) 

XLOM«ABsF(XL0M> 

K0Y*XLAD, and. minus 
IF(K0Y,EO.M1NUS)27,28 
27  CSMlAT«XLAD-XLAM/60|0 
GO  TO  29 

26  COMLAT'XLAD  ♦ XLAM/00’,0 
20  kEYcXLOD.AND. MINUS 

1F(KEY,E0,M1NUS)30,31 
3C  C0Mi.eN«XLe0-XLPM/60|0 
GO  TO  164 

31  C6M10N*XL0D*XL0M/6O|O 
16  4 TlME«RNHR*RNMlN/60 .8 

C ST6RE  JULIAN  DAY,  TIME*  LATITUDE  AND  LONGITUDE  FOR  EACH  Fix 
C 
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1 D ATE  < ISTeflE  >«JUDY 
STImE( IStCRE>=TIME 
STLAT( ISTCRE  >=C0HLAT 
STL0KG(  1ST0RE  )«C0ML0N 
1 ST39E* I ST  SRE*1 
IF (NCwRIT)IOB.  171,106 

r 

r ST£RE  INFORMATION  FOR  WRITING  0N  TAPE 

r 

171  IkRlTE" I WRITE^l 
IYEaR( IwRITE)«IYR 
f'SNTH  ( I I TE  ) * I M 
I5Ay< IWRITE  )*ID 
HR ( JWRI  IE  )sRNHR 
XMlN( IWRITE)«RNMIK*10|0 
2M I N < IWRITE  )*XHIN<Ihr|TE)/10,0 
XLATI  IWRITEJsCGMLAT 
Xl3nG(IwPITE>=C0ML0N 
lFIx(IWRlTE)*iriXTY 
A V F I X < I W q 1 T E ) a I r N 
IF  <2O.IwfiITE>lO9.1O».lO0 

1C  9 MREC*NlREC*20 
\PREC«SPREC*1 

c 

r WRITE  LOGICAL  RECORDS  SN  TAPE 

r 

14  0 -RITE(lOillO)  (CRUISE.  JTMZNE#  IYEAR(  I ) * MONTH!  I ) , 1 DA  Y ( l’j » HR  ( I ) , XM I N ( I 
D.XLATd  ) i XLGNG  ( I > . JFJXU  ),NAVF1X(I  } » I « l IWRITE) 

110  FeRMAT(A8,I5.I2.I2.l2ilXr2,F3,r8.4,F9;4.  !2.7Xl5,24X> 

IF(  IGPT>111, 112,111 

111  IF(  IGPT-2)113, 115,114 

116  F0RMAT  ( XH  ,2XA8(7XlS>,7Xl2.7XI2.7Xl2,7XF2',7XF4.i,6XF8;4‘,5XF9,4,9Xl2 
1,10x15) 

113  IF  (NPREC-D601, 600,®0l 

6 C 0 WRITE (6 1,116) ( CRUISE, l TMZNE, I YEAR (I  ) # MONTH ( I > . X DAY « 1 1 » HR ( I S • ZM 1 N C I 
1)  i XL  AT  ( I ),  XLGNG  ( 1 ) . I F 1 X { I ) ,NAVFIX<  I ) i I»l’,  IWRITE) 

GO  TO  H2 

6C1  IF  (MSD(f.PREC,3)  ,NE.1>  G0  TO  600 
* R I T E ( 6 1 , 5 0 1 ) 
wR  I TE ( 6 1 , 502 ) 

GG  TO  600 

1 1 4 IF (nPREc-1>604, 603,904 

6 C 3 WRITE (6i,  116)  (CRUISE. JTMZNE, l YEAR ( I ) » MONTH ( I ) , I DA Y ( I ) , HR ( I ) , ZM 1 N ( I 
1 ) . XL  AT  C I ) i XLGNG  ( I ) , | F 1 X ( I ) ,NAVFIX(I  )•  I«l‘,  IWRITE) 

GG  TO  115 

6C4  IF  (y,SD(\PREC,3)  ,NE.l)  G0  TO  603 
*RItE(6i,501) 
wRITE(6i,502> 

GG  TO  6 0 3 

115  WRITE (62, 110)  (CRUISE, JTMZNE, I YE  AR ( I ) # MONTH ( I ) , 1 DAY ( i MR ( I ) , Xm!n ( I 
D.XLATd  ) , XLONG  ( I ) , {F } x ( I ) , N A VF  I X ( I ) , I «l',  I WR  I TE  ) 

112  I WR I TE  * 0 
IF(IENDFL)150,108,1>C 

150  ENOFlLE  10 
G0  TO  160 
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r PRecess  the  bathymetry  data  cards 

r 

2 C 0 DC  201  L»1 . 24 
2C1  KE6PHT(L)|0 
KKH  T * 1 
KEEPMT ( 1 ) |MT 
I6ND^L*0 
I wR  l Tc*0 
IFIELD«0 
I81r£C«0 
IBPR6C«0 
i qp  t * i Out 

KNPUTaInPLT 
IF ( I OPT  >241 , 240,241 
IF( 1GPT-2)242.240.242 
WRItE(61,504) 

F5RMAT«iH1,133HSHIP  AND  CRUISE 
1 MINUTE  LATITUDE  L6N&1TUDE 
2ECTED  MATTHEWS) 
wR  I TE ( 6i , 505  ) 

FGRMATdH  ,22HlDENTlFiCATl0N 
1 METERS  ZSNfc) 

D0  202  J*1,ICNTN 
IF(IFIElD)203,205,203 


241 

242 
3C4 


5C5 

240 


TIME  YEAR 

uncorrected 


MONTH  DAY  HOUR 
UNCORRECtED  C®RR 


20NEi65x.44HFATHGMS 


METERS 


read  a Bathymetry  card 


2C3 

--2C4 

2C4 

2 C 6 
2 C 9 
234 

233 


236 


2 C 6 
r2C8 


RE  AD ( 60 , 204  > I SN, JUDY, (BHR(I).BHIN(I)iIDP(I>*I«1»5) 
FORMAT ( 1 2 1 I1,1X,5(F2.F4|1i1XI4) ) 

FORmA T( l2ilXI3,5<lXf2iF8 .1*1X14} ) 

IF ( ISNJ206.207.207 
IF ( ISN*6  )209.208,209 
IF ( I SN*7 >11 , 234,11 
IENDFl»1 

IF  ( 1*RItE)232#232#2<J5 
I0LRFC*lBLREC*lHRlTfe 

I bprec» i bprec«i 

I BEG* 1 WRl T E ♦ 1 

DG  236  K«  l BEG i 20 

I YEAR(K)«0 

mGNTH(K)«0 

ICAY(K)=0 

HR(K)*0 

XMlN(K>sO 

XLAT (K)sO 

XL0nG(K)-O 

NAVF  1 X ( K } 8 0 

ICRMETOOIO 

S l OP  < K > =0 

I F 1 X <K ) =0 

IWRITE>20 

GO  TO  223 

MT* JUDY 

MT«0MR(1) 

CO  210  L«1 #24 

IF ( kEERmT (L ) >211 *21l*  2l2 
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2 1 2 IF  <kEEPmT(U-MT)210»205|210 

21 0 CGNTlNUE 

211  KEEPMT(L)|MT 

GG  T9  205 

2C7  I CNT0«ICNT0*1 
T .ADy*  1 S N • 1 0 ♦ JUOY 
IF IEUD-1 
c 

r COMPUTE  D*Y  AND  month  FrGM  JULIAN  DATE 


CALL  JULlANdYR,  JUDlf.tDilM.LPYR) 
r 

r I\TERF*3LATe  LATITUDE  AND  L BnU J T UOE  AT  THE  TIME  OF  MEASUREMENT 

r 

2C  3 TI“E«BM»(iriEuD)*BMlMlFlkLD)/»0,0 
IF  faMS«C  IFIEL**)  >273,270.273 
27  0 IF(0**I‘.(  Jf  !Eu5)  )273i?7l,?73 

271  IFClC“<iriELC>>273,2C9.273 

273  1 F ( I D aT£ (.  ) • juO Y >202,213, 2i7 

2 3 3 IF(TImE.LE.STImE(j)I?15|?02 

2 1 5 IF(TImE-5TIMF(J))2i7.2i0.2i7 

216  KC«STl9nC(J) 
yC«STl»TI„> 

r ' tfl  t * 

U w | C £ o ' 

217  TNBsSTImE(J) 

TNA.ST  ImE ( J-l ) 

TCX.TJMc 

XNA«STL0NG(J-1) 

YNAsSTLAT(J.l) 

XNB*  STL0N5 ( J ) 

YNB*STLaT(J) 

CALL  PRGPGPITNA.TNBiTDXiXNA.YNA.XNB, YNB.XC.YD) 

• compute  depth  cgrrectibn  accord i no  tb  the  matthews  zone  number 


237  IF(  INPUT. 1)219, 219,216 
210  I DEPTH* I DP  < JFIELD)/1.6208 
GQ  TG  220 

21 9 ICEPTM*IDP( IFIELD) 

220  CALL  MTCBR( l DEPTH , mT . K0RFAT ,K0RMET , METUNC ,MTDC ) 
! WR I T£« I WR I TE*1 


C 

r 


STGRE  INFORMATION  F0R  WRITING  0N  TAPE 


I YE AR ( I wR I TE ) * I YR 
MGNtH ( I wR I TE ) * I M 
I DA Y ( I WR I TE ) * 1 0 
HR( IWRlTE)a0HRC  JFIELD) 

XMlNl  IWrITE ) *BM j n ( lFlbLD)*l0,0 
XLAT(IWrITE)«YD 
XLOnG { l wR I TE ) *XD 
NAVFIX(IWRITE)«ICEPTh*10. 

I CRMET ( IWRJTE)«K0RMfeT 

IF  I X ( I HR  I TE  ) »MT 

IF ( INPUT- 1)274, 274, 27? 


28 
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275  S10P(lWRlTE)«lDPUrJEtD) 

GO  T0  276 

274  SIDP( I WR ! TE  > « IDPt l P 1 ELD ) *1 . 8288 

276  IE (20«JwFJTE)222»222.221 

222  !8LREC«18LPEC*20 
IBPrEC* I0PREC*1 

r 

r hRjTE  LOGICAL  RECORDS  #N  TAPE 

r* 

223  WRITE  1 10.224 ) (CRUISE. I VMZNE, I YEAR <1 ) .MONTH ( I > . 1DAY(  I ) »MR<  l ) .XMlNl  I 
1 > * XL  A T ( I ) « XLCNG  ( 1 ) , NA VP | X ( I ) , ICRMET(  I IFIX<  I l»l«  IUR1TE) 

224  P0RmAT(  A8, 1^,  12,12,  l2ilXr2,r3i?8,4,r9;4,i0Xl5, 15‘,  l5,i6X) 

IP ( I 0 P T >225,226,225 

6 C 7 WRlTE(6l,228  )(CRUlSfc.lTH2Nt»IYEAR(I)#R8NtH(I),IDAY(I>*MR(I),ZMlN(I 
l)iXLATCl),XLONG(I),ZAVrjX(I),SIDP(I)ilCRMET(l)',  IFlX(f )’,  111,  l WRITE) 
22  8 FgRMATCiN  ,2XA8  ,5X I > ,5X  J 2 |5X  12 ,5X  12 , 5XP2',5XF4  ,i  ,4XF8;4-,3XF9 ,4 ,5XF6 
1.1.8XP6.1, 7X15. 9X12) 

GO  TO  226 

225  CO  230  1*1 ,20 

ZK I N C I > =XM  IN  < n/IC.O 

230  Z A V F l X ( I ) ■ \ A V r IX(I )/l0, 

IP ( lGPT-2)227, 229,231 

227  IP(  I9PREC.1)605,607I605 
6 C 5 IP <MSD<  I9PRSC.3) .NE,1)  GO  T0  637 
wR 1 TE ( 6 1 , 5 04  ) 

«R l T E ( 6 1 1 5 05  ) 

GO  TO  6g7 

231  IP ( I0PPEC*1)6O9,6Q8|6O9 

6 C 9 IP(mCD( IBPREC.3) ,NE|1)  Go  T8  608 
wR  I TE { 6l , 504 ) 
kR I TE ( 6l , 5 05 ) 

GO  TO  608  .... 

6 C 8 kRlTE(6l,228  ) (CRUISE. ITMZNC, iyEARI l ) * MONTH ( I ) , I DAY ( I ) , MR ( I ) , 2M 1 N ( I 
D.XLATd  I.XLONGdl.iAVFjxdJ.SIUPdJ.tCRMETd),  If  IX  < ! > , I «1 , l WR  I TE  > 
229  «RItE(62,224  ) (CfiUlSt. JTMZNE. I YEAR (I ) .MONTH d > , I DAY ( l ) . HR d > , XMI N d 

l).XLAT(I).XL0NG(I).NAVrix(I  ).ICRMET<I)',  iflXCl)'.  1 «1  • 1 WRt  TE  > 

226  I HR  I TE*0 

IF  < IENDTL  )232 , 221,232 

232  ENDPlLE  10 
GO  TO  160 

221  IPIELD*IPIELD*1 

IP ( IPIElD«6  )203,205|205 
2C2  CONTINUE 

mR  I TE ( 6 1 ,233 ) . . . . „ 

233  FORMAT ( lHO  .44HUNA8L'  T®  INTERPOLATE  LATITUDE  AND  LONGITUDE) 

STOP 

r. 

C PROCESS  THE  MAGNETICS  SlATA  CARDS 

r 

3 C 0 M * 1 

lOPT^JOuT 

IENDFL'0 

IHRITE-O 

M'.RECpO 

MPRECuO 

IF  SELD-0 
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I F ( I OPT  >341,34  0,341 

341  IF ( I0PT-2)342,34O,342 

342  -RlTE(6i,507) 

5C7  F6RMAT<iHi,i29H$HI*  aN0  chuJSE 
1HSUR  MINUTE  LATITUDE 
2U*l  MAGNETIC) 
wRJ’E (6l»5O0) 

5c8  FCRMAKiH  ,23HlDENT|FjCATleN 
INTENSITY) 

340  D0  301  J*1 i I CNTN 

IF(IFIElD)303,302,303 

p 

C READ  A MAGNETICS  CARD 


time  year  month 
LONGITUDE  TOTAL  MAGNETIC 


Z0NE',73x,29HlNTENSltY 


DAY 

R8SID 

I 


Z3C2  READ<06,304)JUOY,(RMHR{|),RMMIN(I),RMSEc(I).ULTml!MO(!),l 

X i al 1 12 ) 

"3C4  FeRMAT(I3,6(3F2,I4,l5)/8{3F2.l4,I5)) 
r-  IF  (EOF, 06)305, 306 

C3C5  IFLIGhTsIFLIGhT 
' IF ( IFlIGMT-13)302,302#325 

3C2  READ(60,304)ISN,JUDY.RMMR  . (IMOCI  ),Ui;i2i 
3C4  F6RMAT(I2,I3,iXF2,lia)(|5)) 

IF ( I Sn  ) 305 , 306  * 30  6 
3C5  IF ( IS N *8)11, 325, 11 

325  IF.NDFl«1 

IF  < IWRITE)324,324,326 

326  MLREC*MLREC*IWRITE 
MPREC^MPREC*1 

I BEG*  1 WRI  TE*1 

L0  327  K«  I BEG ,20 

I YEAR  (K  ) ■ 0 

MONTH  ( K ) ■ 0 

ICAy(K)*0 

XLAT(K)xO 

XLOnG(K)«0 

MR ( K ) « 0 

XMlN(K)=0 

ZM I N ( K ) s 0 

NAVFIX(K)«0 

327  ICRmET(K)»0 
I HR  I TE*20 
GO  TO  317 

3C6  ICTMAG*ICTMAG+1 

ifield«i 

C 

£ COMPUTE  DAY  AND  MONTH  FROM  jULUN  DATE 


CALL  JULlANUVR.JUDY.lDjlM.LPYR) 


INTERPOLATE  LATITUDE  and  LONGITUDE  AT  THE  TIME  OF  MEASUREMENT 


303  TM«ITM( (FIELD) 

TlME«RMHR*TM/60 ,0 

0303  TM*RMM I N ( I F I ELD ) * RmSEC ( J r I ELD  )/6C , 0 
c T I ME*RMHR  < I F I ELD ) ♦ TM/*OiO 

C IF( IFlIGHT-2)7O0.7O0,70; 
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1 

1 

1 


C7C0  ALT«1200, 

7 G0  T®  702 

rlil  ALTi!ALT<  1 r I ELD } *100 

If  UHG<  1 FIELD) ) 350. 3l>  1^50 
350  lF<lDATE<g>-JUDY)30l.J07.3oe 

3 C 7 1F(T!ME.LE,ST1ME(J))309|301 

3 C 9 IF (TlME-STlME(J})30»,3l0.3oa 
3 1 P XD«STL0NG(J) 

YD«STl*T(w) 

G0  TO  330 
3 C 6 TNBsSMmE»J) 

TNAsSYImE(J*1) 

TCX.TJ.'E 
XNAsStLGN'G  ( J-l ) 

YNA»STLAT(J*1) 

XNB»STL0NG<J) 

YNBsSTLaT(J) 

CAUL  PR0PER(TNA,TNB|TDX|XNA,YNA,XNB,YNB,XD.YD) 

CALCULATE  RESIDUAL  MAGNETIC  INTENSITY 

CALL  SPhERE(xD.YO,GLSNiGlAT) 

ZYR«IyR 

IF (lPYR-2)311, 329,311 
LEAP«366 
GO  TO  3l2 
LEAP«365 

LATE«1900,  ♦ ZYR  ♦ JFL0ATF(JUDY)/FL0ATr(LEAP)) 
ITYPEsl 

IF ( ALT ) 31 A 1 313  # 314 
ALTsl.O 

BLT*,30a6  • ALT/1006.  . , 

CALL  IGrF(DATc,ITYPE.ULT.GLAT,OuGN,XV,YV.ZV,TV) 

I T V « T V 

:OP*lMG(IFIELD)*ITV 
STAlT|ALT 

C ST6RE  INFORMATION  FOR  WRITING  0N  TAPE 
C 

1 WR I TE ■ I WR 1 TE*1 

1 YEAR ( I wRJ  TE ) si YR 

M0NTH(  IwRITE  )»IM 

ICAY(IWRITE)=ID 

XL AT ( I KR I TE  ) * YD 

XL0NG( IwRITE)*XD 

HR ( l NR  l TE  ) sRPnR 

XMlN( IWrITE>sITM< IFIELD)*10 

ZMINUWRITE)*ITM<IF|ELDI 

NAVFiXt  1WR1TE)*IMG(IFIELD) 

1CRmET(  ILRITE)«NDP 
IF (20. IwRITE>316, 319,31? 

3 1 6 MlREC?MlREC*20 
mPREC5HPREC*1 
r 

c NR  I TE  LOGICAL  RECORDS  0N  TAPE 


C 

330 


329 

311 

3 1 2 


313 

314 
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31? 

3 1 • 

31* 

321 

611 


322 
6l  3 


610 

612 

321 

323 

320 

324 

319 

*319 

3C1 


KH|TK(iO«9il)(Citt|tk»|Tllfll6*tVtAW(t)«NMtKlI>f tOiVMVa  hR(I>.xmjn( 
11)»xLaT(  I > , XL  CNG  < X )|Mvr  iXcn.lCRXETd),  j «i . | wa  1 TE  ) 

IF (iepT)Sl9.320.3l« 

If ( I OPT -2 ) 32 i *323,32? 

IF(hRREc*1)6;1.610.6U 
|F(*5D<MPR5C.3>.NE.i>  0*  *3  610 
WRlTE(6l|907) 

«R  I TE ( 6 1 •90S  ) 

Gc  TS  6 i0 

I F ( KPREC*1 ) 613,612(613 
IF (M0D(mPREC,3) ,NE.l)  00  T 3 612 
WR  I TE ( 6i*9  07  ) 

-RITE (61,930) 

G0  TO  6 1 2 

«R,TE<01,320  )(cRUISS.JTHzNE(,YEA<»(I),MeNfH(l).IOlY(n#  HR ( I ) • ZM l N ( 
li ).XLaT(I  ).Xc£NG(I >|KAVF|Xcl > . I CRHET ( | I . f « 1 • I u« 1 T£ * 

G6  T t 320 

-RlTE(0ij320)(CRUISk.,TH7NE,IYEAR(I),M0NfMCI>,!OAY(r»<  hR(!).|min( 

XiilX^AT(  1 ) *XL«NO€|  > «6<AvriXM  > . IC^HETC  1 9 . j • 1 . | W« I TE  > 

'2x*ft‘6*,>'’>f*5,|2,'**r2'®*r«.1.5*F0.4.«xr*:4,,ax!9. 15X15) 
wRI TE (62 ,310  MCRL  JSt . I TMjNE , I YEAR ( I ) , MOnT H ( I 1 . 1 DA Y ( I $ , HR( I } ,XMlN< 
II 1 iXLATc I ) ,X10NG( I ) , KA Vf I X C I )*  ICRHETt I ),  ill, I WRITE) 

U*I  TE«0 

IF  MENDFU324.3l5.324 
ENDFIlE  10 

gc  re  160 

I F | ELO« I r I ELD  ♦ 1 

I F I ELD* I F I EL0* 1 

IF ( IF  JEl 0*13)303 #302, 303 

CONTINUE 

wR I TE ( 6i • 233 ) 

STOP 


WRITE  0UT  ACCUMULATED  DATA 

4 C 0 WR I TE  ( 6 1 , 5 1 0 ) 

910  FGRMAT(IHB) 

WRJTE{6i,401)CRUISE 

401  F£RmAT 'i^1.35HSHIP  AND  CRUISE  I DENT  I F I CAT  I ON  ,A0) 

IF(NGWRIT)  173, 172,173 

172  wR I TE ( 6i • 403 ) 

4C3  F£R«AT(iho,15HNAVIGATJ(5N  DATA) 

WRITE(6i,404)ICNTN 

404  F0RMATJ1HO.34HNUMBER  OF  NAVIGATION  CARDS  READ  ■ .14) 

WRUEtei,  405  >NLREC 

^6RmAT(iHo,44mNUMBEK  OF  LOGICAL  RECORDS  WRITTEN  ON  TAPE  ■ ,14///) 

173  IF  (N0wR(T«l)i74,174|l79  * 

174  WR I TE ( 6 1 , 407  ) 

4 C ? F5RmAT(iH0,15HBATHYMFThT  DATA) 

IF(KNPUT«1)408,408,4o* 

4C8  WRlTE(6l,410) 

410  FORMAT(iHO,39HDEPTH  DATA  GIVEN  IN  UNC0RRECTED  FATHOMS) 

GO  TO  4 il 
40*  WR I TE ( 6l , 412 ) 

4 1 2 FORMAT (1HO,30HDEPTH  DATA  GIVEN  IN  UNC0RRPCTED  METERS^ 
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411  WRITE (6 1,413)  (KEEPmT(R)|K«1,KKHT) 

413  r£RM4T<lH0,35HH4TTMtwS  20NES  PASSED  THROUGH  WERE  ,J4f4) 
wR|tE(6i,414)  ICNT0 

41 4 F£RmAT<iHQ,34HNUm8EH  tir  BATHYMETRY  CARDS  READ  • ,10) 
wRItE(6i,415) I BlREC 

4 1 5 PeR«AT(lW0,36h\U^9E«  BE  lEGICAu  RECGRDS  WRITTEN*  ■ 1 6 / / / J 

175  IF EwRlT.l  ) 176, 177,176 

176  wRITE(6i,417) 

4 1 7 FCRm A T (iM0,14hm4GNETICS  DATA) 

»RlTE(6i,410)STALT 

410  r£«MAT(iHo,43HW*E!GHT  IN  FEET  A90VE  0R  BELGW  SEA  LEVEL’  IS  ,F3> 
hRlTE(0i,4l9)ICTMAG 

4 1 9 f£RMAT(iHo,33HNUM3ER  0F  MAGNETICS  CARDS  READ  » ,ll?> 
wR  I TE ( 01 , 4 20  )MLREC 

420  F£RmAT<iH0,36MNUMBER  6F  LOGICAL  REC0RDS  WRITTEN  ■ 1 1 2 ) 

177  ENDElUE  10 
REWIND  10 
STOP 

END 


i 

1 
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5.4DS  Gc  £ TA  T A 

IDEMT  GEODATA 

sRfGRA-  LENGTH  J3*05 

cNTfiv  PGlNTS  GEGDATA  do 76 7 

EXTERNAL  SYMBOLS 

C 5 GENTRY 
TREND, 

C80STGPS 
Q1C1C100 
C 8 0 LM  C T , 

SKIPFILE 

backfile 

JULIAN 

PRGPQH 

MTC£R 

SPHERE 

!GRF 

xkolf 

cao i fest 

EFT, 

REw, 

TSh  , 

STH, 

CNSINGL. 

00625  SYMBOLS 
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SCPRGjTUE  JULIAN  ( I Y'»  J jnY  . 1 0. 1 M . LPY«  > 

- C\VE*tS  JULIAN  Sate  I NT0  Jay,  month,  year  for  ANNOTATION  IN  3CEAN0 

' IT  » YEAR,  JtCY  « JULIAN  Day,  ID  * CALENDER  Dav,  1M  ■ CALENDER  MON 

- cpplIments  of  soe  FtrtN  - code  e 1 7 4 • 17  nov  1571 

c ^EDIFIED  by  LE2N  LA  LUMjF«E  - CODE  8178  • 17  NSV  197? 

r 

Cl^ENSlGN  JJ(13),  A A ( 12 ) 

type  integer  aa 

CATA  (tjJCI).  1 « l!  13)  » o.  31,  59,  9 o',  120,  l5l',  lSl,  2l2,  243, 
1273,  304 i 334,  365) 

DATA  ( ( a A ( I ) , I * i,  12)  • 1,  2,  3,  47  5,  A,  7,  8,  9.  io',  11.  12) 

LL  • 1 Y 
rf* JUDY 

IF  (LL.LT.O)  GO  T s 10 
IF  (LL.EC.O)  GO  T0  1 
IF  (H0D(LL,4)  ,EQ.O>  GU  TO  4 

1 IF  (MM.lE, O.OR.MM.GT  ,362)  G0  Tfi  4 

2 CO  3 K>2 , 13 

IF  (Mm.lE, JJ(K)  ) GO  T«  6 

3 CONTINUE 

4 uFVR»2 

IF  (MH,lE,0.OR.KM,GT.366)  GO  Tfl  10 
IF  I M“  , lE  , 59  > GO  TO  2 
IF  (MM,EC,60)  GO  TO  5 
**  “ a M « 1 
GO  T 3 2 

5 1 1 -2® 

K *3 

GO  TG  7 

6 I I*uM»Jj(K„l) 

IF  (LL.EQ.O)  GO  T0  B 

IF  (P0D(LL,4)  ,NE.O>  GO  TO  8 

IF  (H0D(LL,4>  ,EQ.0.AnD.Mm,lE.59)  GO  T0  0 

7 IF  (LL.GE.10)  G0  T9  8 

8 I C» 1 1 
if-sAA(K-l) 

6 RETURN 

10  PRINT  11,  LL • MM 
GO  TO  9 

r 

11  FORMAT  (1X,*YEAR  ■ •.  l2ilX,*JUL!AN  D*Y  » •'<  1 3 , ix',  • ARE  BAD  OATES*  ) 
END 


F 1 

r 2 

F 3 
F 4 
F 5 
F 6 
F 7 
F 8 
F 9 
F 10 
F 11 
F 12 
F 13 
F 14 
F 15 
F 16 
F 17 
F 18 
F 19 
F 20 
F 21 
F 22 
F 23 
F 24 
F 25 
F 26 
F 27 
F 28 
F 29 
F 10 
F 31 
F 32 
F 33 
F 34 
F 35 
F 36 
F 37 
F 38 
F 39 
F 40 
F 41 
F 42 
F 43. 
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9.4DS  .ILIAN 


IDENT 

OQAM  LENGTH  00301 

CNTflY  PSlNTS  JULIAN  00053 
PXTEPNAL  SYMBOLS 

THEiSD , 

CBCuICT , 

XM0DP 

STH. 

ONSINGL, 

00116  SYMBOLS 
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SlBRSuTJNE  PR0P0R  (TnA,TnB#TDX.XNA*YNA,XNB»YNB,  XD»VD*J  H 1 

r EPPUTES  THE  X-Y  C80HrjNATfcS  F0R  PLATTING  BATHYMETRY  M 2 

r PROGRAMMER  - LE0N  la  LURlfcRE  - CADE  «1?4  • 15  DEC  1971  H 3 

- H 4 

pERISD«tnb-Tna  h 5 

CELT«T0X*TNA  H A 

RATIS«OELT/PERI0O  h t 

CELX«XNB»XNA  H 8 

CELY«YNB»YNA  H 9 

xC«DElX.RATIG*XNA  H 10 

YD«DELT*RATle*YNA  H 11 

RETURN  h 32 

END  « 33. 
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9.4DS  PfioPSR 


length 

cNTRy  POINTS 

EXTERNAL  symbols 


PRGP0R 


00C54 


SYMBOLS 


OflODlCT, 


IDENT 

00140 

oooo  j 


PROP0R 
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FTN5.4A 


SUERSUT INE  SPHERE (CIGN.OL AT. GLQN.GL AT ) 

IF  <010^)10,20,20 
GL0N«DL0N 
GG  TO  30 

Gl0N*36 0 • ♦ O10N 
Gt AT*90  . «DL* T 
RETURN 
END 
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SUBROUTINE  MTC0R  < I l)EP . MT  ,KCDF  , KCDM , METUNC',  MTDe  > 

JAN  10,  1972 

MED  I r I ED  F0R  CDC  3800  BY  LEON  LA  LUM l ERE  • CODE  8174  • 20  MAR  1972 

SUBROUTINE  MTCPR(AP(i)i  MATTHEWS  TABLE  CORRECT  1 BN 

VERSION  WHERE  COtFf  lC  JENTS  entered  AT  RUN  T J MB 
alculates  corrected  depth  in  meters  when  given 
MATTHEWS  Table  coefficients  and  uncorrectcd 
OEpTh  IN  FATHOMS', 

COPIED  FROM  WHO  I PR9GKAM  FOR  IBM  1130 

input, 

IDEP  « UNC0RRECTED  DEPTH  IN  FATH0MS 
MT  ■ MATTHEWS  Table  CflEFFJCIENTll  • 52 1 

OUTPUT, 

KCdF  ■ CORRECTED  DEPTH  IN  FATHOMS 

KCdm  ■ corrected  DEPTH  IN  METERS 

MTdC  » MATHEWS  TARtE  DEPTH  C0RRECTJ8N  (METERS) 


internal 

nudf 

wudf 


NUDF  WATER  depth  U::CORRECTED  fathoms 
wudf  WATER  DEPTH  UNtORRECTED  FATHOMS 
WUDM  WATER  DEPTH  UNtORRECTED  METERS 
wC0fl  CORRECT  I 0N  IN  METEHS 


SETTING  MATTHEWS  TA0LE  COEFFICIENTS  SUPPLIED  BY  C,  GANTAR 
DIMENSION  mTn(52),AMT(52)iBMT(52),CmT(52),DMT(52),EMT(52>,FMT(52) 

REAL  IDEP 

DATA  (MTN. 1,2, 3. 4, 5, 6, 7 ,8, 9, 10, 11, 12, 13, 14, 15,16,17, 16^19, 20. 21, 22 
1. 23, 24, 25 ,26. 27 ,20, 29> 30.3i.32, 33, 34, 35, 36, 37. 38. 39, 40. 41,42.43,44 
2.45,46,47,46. 49,5o, 9n52)  , . 

DATA  (A mT*-,05. ,21,11*,??*,, 35, ,43. .12, -.21, .76, *,19, .,o7,1.29, ,34 
I'*i04,,l2,.39,i,i6,1.25j-|ll,,44,-,31r,  .4  3,-.29,-.06,-,110i»,08,”.2 

23. . 06. 1. 14. 38.. 2. 34,  .59, 1,15,  ,58',  .01,  .59, 2^49, .97,3,93.1.9 
36,*.15,.09,I06,I42,*.05|.03,-I4i,-,19;-Ioi',  .07) 

DATA  ( BhT»-,0  1683, 00247,  . 006885,  , 000632.  . 015719,  , 004269,  .013881, 

1. 025565. . 0307 d9, .023035, . 03 1673,. 03649 j, .031972,, 049784,. 03225S, . 0 
24  0042,  , 037533.  . 0435*1  ,.  Q46531,  ,034097; , 0i9659, . 025f|9S.  , OJ.2063,  . 0 05 
3933  ,*  , 000139  . - . 0C637<> -j  COR94,-.  012026', -*01*3l7,-,  oil  231  , . 03*103  . . 
4 037121,  .03787  4.  .03524-,  . 031429,.  0378  72,.  0*3038,  . 03*5 , . 0257  75 , . 03583 

53. . 02 1001. .0186. .0219  06, .013118  , .001109,. ' 0*4699, ,02651, .02*0 06  , .03 

•5524, ,038*1,. 049807, .000832) 

DATA  ( CmT s.  1 3 0*95  , -,c2fl 055, -.  033 4 16  , ,0307  85  ... 169795, ,0494 19, -.059 

1931..  , 156579,  - .1412B6-,  . 00*79  10  1806',  l 2 1 8 7 1 ',  -,  0*0565  ,..  233986  , . 

2 011260  , - ,363995  , -.227630...  239633  ,-.383832'. -.2*561', -,0*4804, -.2165 
3o9,-,014753,-.0  0920*.io2l74  6. . 0 82109, Ti oi 101. .118*5,.  193882, .0632  7 
48, -. 255508, -. I747i2»-i 1*6*92. -.149662,-. 1*1422, -,178i2,-. 210*77,-. 
5103*53, -.0*2957, -.1*2494,-.  0*7*1, -,o227i', -.117073,-,  S39335,  .023648 
6, -,  27 l7o5,-. 0*1812, i 0**156, -.059269  ,-; 097981  , . 001276t, 161859) 

DATA  (DhT*-.  0*0587,  , 017161  ,-.  00*71  ,.  012525’,  . 1 1527  i',  -,  015992,  . 0*609 
14,  , 083656  , . 057  27  ,..  0i9  41.  , 036137  , , 03*7  75’,  T 00*688 , , o83*87',  -.  030352 , 
2.207  727  , ,1027  7 4 , .0*91 35,  .187189,  .123539,  . 022*5'.  . 119*29', . 00»ll8  , . 02 


J 10 
J 11 
J 12 
J 13 
J 1* 

J 15 
J 16 
J 17 
J 18 
J 19 
J 20 
J 21 
J 22 
J 23 
J 2* 
J 25 
J 26 
J 27 
J 28 
J 29 
J 30 
J 31 

J 32 
J 33 
J 3* 
J 35 
J 36 
J 37 
J 38 
J 39 
J *0 
J *1 
J *2 
J *3 
J ** 
J *5 
J 46 
J *7 
J 48 
J *9 
J 50 
J 51 
J 52 
J 53 
J 5* 
J 55 
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2 


r 


r* 


3^232,  .Olfill,-.  013315. -.020412, -.030545',-;  061797. -,0  01626,.  121471,  . 

4064 336. . 081001. .056314,  j 0^5191, ,067, ,083684  , , 020131 , ] 037603 , .092  35 

59. . 018840.. 0092 4. .000 13 1.023 7 95. .007413.. 10 30 95  , ,167905,, 004016, .o 

670316..  097068, -.000464, *.162070) 

DATA  (tMTp. 092394,, 04bO95(. 127138, •,02266»». 24916  6',  .0  69  75, ..007349 
1* -,154994,-, 075775, ,oBl0o4,«. 035675, -704Q077, ,023167V* , It 59 05, ,io7 
256  , 47  0139  , -,  166568,-,  137270,..  366057',..  235691  ...  029841, -;232503, 
3. G15910, ■,046104. «.Q4S777i,  02  4580, ,o37o58». 002062#, 110005 1 '.001406, 
4-, 2l49»2i«. 00685, -,13Bl37|.,o79i56,-, 07510',  - ,09377,  - , 136212,  .01265 
52 0 41351,-. 106534, -J014159, ,001108. -To9774,-, 025665. -.004717, -.1 
64274  0 ,..903392  , . 000155,  *.170839, -.2527',  .24145,  ,983392) 

CiTA  <^74-, 072171, *.150140,  -.231021,012058,,  190  970,.  ,070755..  004  6 

179..  109779,.  030023,  *.07  5341,. 011294,  , 015287 ,..  02  7333*.  , 062495 ...  096 

2216,  , 403122,  .100  73  , i 074903,  ,267798  , . 1692  77',  . 021495',  ,16502. -.  028903 

3. . 0  35913, ,0  40604,., 0l4l77.., 022102, -,045831... p6729, 7001066,  .1410  0 
47,  , 04960,  ,094  743,  ,04808  5.  ,039534  , , 0 50  990',  VO 89176',. ,032552,  .017001  , 

5.140729  ,,  003067 , - , 0 04373 , 5 050355  , , 0 1 05l5',  7000102 ',  , 067775', 2 . 003205 , 
0-, 00155.  ,149654  , .235294 1-!623200  , -2 .003205  ) 

NZERQtO 


NUDF  * I D£P 
I F (NUDF)  3,5,1 
1 IF  (52-MT)  4.2,2 

WODF-IDEP 
wIDm*wUqF  *1,0208 

wC0R«AMT(f»T).BMT(MTj*liUOM*CMT<HT>*l.E-D4i{HUDM4*2)*OMT(Mt)«l.E-O7* 
l(WUDH..3).EMT(MT)*l,Efll.(WUDM**4)*FMT(MT>*l,E.l5*([WljDM**5) 
^CDM*WUDM.UCQR 
KCDm*wCdm-0 .5 
f'TDC*wC3R*0.5 
KCDF*(WCDM*O.5408l)*Of5 
METu^C*WUCH 
RETURN 

3  PRINT  0 
G0  T 0 5 


4 PRINT  7 

ERROR  Exit 

5 kCDm*NZERS 
MTDC*NZER6 
KCDF  *NZER0 
RETURN 

6 FORMAT  (IX, ‘DEPTH  READ  IS  NEGATIVE*) 

7 FORMAT  (1X,*MATTHEW6  ZONE  NUMBER  IS  GREATER  THAN  52*\ 
END 


J 56 
J 57 
J 58 
J 59 
J 60 
J 61 
J 62 
J 63 
J 64 
J 65 
J 66 
J 67 
J 68 
J 69 
J 70 
J 71 
J 72 
J 73 
J 74 
J 75 
J 76 
J 77 
J 78 
J 7** 

J 80 

J 62 
J 83 
J 64 
J 05 
J 66 
J 87 
J 08 
J 69 
J 90 
J 91 
J 92 
J 93 
J 94 
J 95 
J 96 
J 97 
J 98 
J 99 
J 100 
J 101 
J 102 
J 103 
J 104. 


42 


NRL  REPORT  7861 


5.4DS 


PTC~R 

I DENT  MTC0S 

3RJ GRA y LENGTH  0l0lS> 

= MRy  PGInTS  MTCGR  00577 

EXTERNAL  STHB3LS 

01Q10100 

THEND, 

OflQDlCT, 

STm 

00074  SYMBOLS 


I 


mo  i o .1  o n .o  o 


BLODGETT  AND  MASSINGILL 


SUBROUTINE  I GRT  ( date  . ITYPE.ALT.COLATiELONC.X.  y’.Z.T) 

-'GENP  I G«F 

r 

r* 

r this  Subroutine  computes  values  or  x.v.z  and  t for  a given  epoch 

" ANC  POSITION  FRGM  THE  SPHERICAL  HARMONIC  COEFFICIENTS  OF  tH8 

o International  geomagnetic  Reference  field 

C Ca7a  DATE  ■ REOUIRED  fcPOCM  In  YEARS  AND  DECIMALS  OF  A YEAR*,  A.O. 

C I TYPE  * 1 IF  GEODETIC  COORDINATES  are  BEING  USED 

r I TYPE  ■ 2 IF  GEOCtfjTRjc  COORDINATES  ARE  BEING  USFD 

C ALT  ■ HEIGHT  ABOVE  MEAN  SEA  LEVEL  In  KIlOMFTERS 

" AlT  a RADIAL  DISTANCE  FROM  CENTRE  OF  EARTH  IN  KILOMETERS 

c Col at  • colatitude  in  degrees  (o.o  to  leo.O) 

r ElONG  ■ EAST  LONGITUDE  IN  DEGREES  (o.o  TO  360.0) 

r,  OUTPUT  X * NORTH  COMPONENT  OF  MAGNETIC  FORCE  IN  GAMMAS 
Y » EAST  COMPONENT  PF  MAGNETIC  FqrcE  IN  GAMMAS 
2 - VERTICAL  COMPONENT  OF  MAGNETIC  FORCE  IN  GAMMAS 

(POSITIVE  DOWNWARDS) 

T * TOTAL  MAGNETIC  FORCE  IN  GAMMAS 
N.B.  The  COORDINATE  SYSTEM  FOR  X,Y  AND  Z is  THE  SAME  AS  THAt 

specified  by  itype 

DIMENSION  P(44),G(44),cL(8),SL<8), AGH(flO) ,DGH(80) 

The  S.h.COEFF IEIENTS  AHe  SUPPLIED  BY  The  FOLLOWING  DATA  STATEMENTS 
ALTERNATIVELY,  TmE  COEFICIEnTS  MAY  BE  SUPPLiED  BY  STATEMENTS  OF 
THE  FORM  - AGH(i)  * .30339*  AGH(2>  * -2i23^  ETC 
OR  THE  COEFFICIENTS  MAY  BE  RkAD  IN  THE  MASTER  PROGRAM  AND  MADE 
AVAILABLE  TO  THE  subroutine  By  A COMMON  STATEMENT 

e* 

0 SET  INITIAL  VALUES 

r 

data  (Agh  3 -30339  . 2123  , , 5758., -1654  ;, 2994, ,.2P06,’. 1567. ,130  ., 
Al  297  <# -2036.  ,-40  3 . , 12891 ,242.  , 843,, -176  .,  958,  ,S05,',  149,  , 492  ,, 
B-280|l.392.  , 8.  , 25  6 , ,-26  5.  ,-223,,357.  , 16.  , 24  6 , , i25 , ', -26  , . 123 . , 
C-l^i, , • 107 , i -5l , ,77  , , 4 ^ * . P 0 , , • 1 4 . 1 4 , , 106, , • 2 2 9 , #68, , 3 , #-32.  , • 4 . , 
C-l0.i*il2,,-l3,,7i.|-54if,»57.,0.,-27,,i2,,  -8, ,-25, #9, ,-9. ,23,, 
El3,,*l9,,»2.,-l7.ll0.i9f.3.l'3,,-l3,»-i2.>5.#-4.,«i7,,7,,4,,»5., 

F 22 1 ,12, , *3 , #6, , -16, ) 

DATA  (CGH  « i5.3,8,7,*2j3,-24,4,o,3,«li.8i»lI6,*l6,77o,2,*l0.8, 
A4,2,0,7,0,7|-3.8,«7|7',  •C.7,0.2,-O.J#*3.0,1.6»”0,li2«9,»2,l»-4,2, 
Bl, 9, 1,1, 2, 3, 2, 9, 1.7, o',  6, -2. 4, 0,0,0, 8,113'.  - 0.3,-0,l',-8,3,»0;9,l.l, 
C*0,4(l,9,2.0,»0,4l-l.l,«o,4,o.l,-0,2i0.9.-0.5,«0,3,»^,l,»0.7,0.3, 
C-0,5,0,4#o,3,0.2,0.y.0.4.»0.2l0.2,»0,6’,  0 1 3*.  0 . 1 ',  0,4',  0 , 1 0.6, -0.2, 
E0,O,»0i3i0.O,'O,2,<'0.1,»Oi3iO#3i-O,4,tO.3#*Oi3i-O,5,«0,3) 
T«DATE-i965, 0 
R*AlT 

ONErCOLAT *0,0174533 
SLATfCOs(ONE  ) 

ClAT»S!N<0\E  ) 

0NEsEl0nG*0, 0174533 
CL ( 1 ) »CqS ( ONE  ) 

SL ( 1 ) sS I N (ONE  ) 

XsO.O 

ysO.O 

z«o.o 

C D ■ 1 , 0 
SC«0,0 


JS170700 


JS170800 

JS170900 

JS171000 

JS171100 

JS171200 

JS171300 

JS171400 

JS171500 

JS171600 

JS171700 

JS171800 

JS171900 

JS172000 

JS172100 

JS172200 

JS172300 

JS172400 

JS172500 

JS172600 

JS172700 

JS172800 

JS172900 

JS173000 

JS173100 

JS173200 

JS173300 

JS173400 

JS173600 

JS173700 

JS173800 

JS173900 

JS174000 

JS174100 

JS174200 

JS174300 

JS174400 

JS174500 

JS174600 

JS174700 

JS177900 

JS178000 

JS178100 

JS1702OO 

JS178300 

JS178400 

JS173500 

JS17S600 

JS178700 

JS178800 

JS178900 

JS179000 

JS179100 


44 


NRL  REPORT  7861 


L«1 
Hit 
N • 0 

G0  T9  (1,2),  I T YPE 

r 

r VERS  1 0N  FROM  GEGDETjC  T0  GEOCENTRIC  COORDINATES 
C 

1 A2M0680925. 

32*40408585, 

0NEsA2.ClAT.ClAT 
TW0sB2«SLAT*SUAT 
ThREEs3\E*TWS 
FSUR*SQrT (THRfcE  ) 

«»S:R T( ALT •( At T42,o»F0UR)^(A2*flNE*B2*TW0J /THREE! 

CC»( ALT*FSLR)/R 

SC«(A2-B2)/F0uR»SLAT«CL*T/R 

3NEsSlAT 

SLAT«SLaT*CD-CLAT*SC 

Cl*T*CLaT»CD*0NE«SD 

r 

2 RATIO. 6371.2/R 

r 

£ COMPUTATION  0F  SCHMIDT  0UASI-N0RMAL  COEFFiclENTS  P AND  X(.Q) 

P ( 1 )*2 ( 0»$LAT 
P ( 2 ) *2,  O’CtAT 
P(3)»A,5.SLAT.SLAT-i.> 

P(4)»5,1961524.CLAT*SLAT 
Q ( 1 ) • «Ct AT 
0 ( 2 ) *SLaT 

Q(3)»-3.0«CLAT«SLAT 

Q(<)«1,7320508*(SLAT.SLAT-CLAT.CLATJ 

C0  15  Krl,44 
IF(N.M)  3,4,4 

3 H . 0 
N*N*1 

RR.RAT 1 0. • (N*2 ) 

FN*N 

4 FM«M 

I F C K • 5 ) 0,5,5 

5 IF(m.N)  7,6,7 

6 0NE;SQRT(1,O-O.5/FMJ 

P(K)«(1.0*1,0/FM).ONE*CI,aT.P(  J) 

C(K)B0ME»(CLAT*Q( J)*SLAMm.P( J1 ) 

SL(M)sSl(M-1)*CL(1) *CL (M-l )*SL( 1 ) 
Cl(M).Cl(H-1).CL(1)«SL(M-X)«SL(1) 

G0  T0  0 

7 ONEsSQRt(FN.FN-FM.FM) 

TH0sSQRT( (FN-1 . o ) ••2-Fh»fM)/0NE 
THREE»(2,O.FN-1,O)/0nE 

I *K-N 
J*K-2.N*1 

P ( K ) F ( Fn*1 , 0 ) • ( THREE »5L At/FN.P ( I )-TH0/(FN»i , 0 ) *P ( J) ) 
O(K).THaEE*(SLAT.Q(J).CtAT/FN.P(I))»TW0.O(J) 


JS179200 
JS179300 
JS1794  90 
JS179500 
JS179600 
JS179700 
JS179800 
JS179900 
JS180000 

jsieoi oo 

JS180200 
JS160300 
JS180400 
JS180500 
JS180600 
JS180700 
JS1808  00 
JS180900 
JS16100  0 

Jsienoo 

JS181200 

JS181300 

JS181400 

JS181300 

JS181600 

JS181700 

JS181800 

JS181900 

JS182000 

JS182100 

JS182200 

JS182300 

JS132400 

JS182500 

JS182600 

JS182700 

JS182800 

JS182900 

JS183000 

JS183100 

JS183200 

JS183300 

JS183400 

JS183500 

JS183600 

JS183700 

JS183800 

JS183900 

JS184000 

JS184100 

JS1042OO 

JS184300 

JS164400 

JS184500 

JS164600 

JS184700 
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SYNTHESIS  OF  X,Y  AND  Z IN  5E0CENTR1C  COORDINATES 

ft  GNEs(AGh<L>*DGH(L>*T)*RR 
I F < M ) 10*9.10 
9 X«X*0NE*O(K) 

ZaZ-ENE*P(K) 

L*L*1 
GG  T9  14 

10  ThS:(aGh(L*1)«DGH(L*1>*T)*RR 
THREE s0\E*CL(H)*ThG*SL<M) 

XsX*THHEE*Q<K) 

Z«Z-ThR££*R (K  ) 

IF  < CL AT  ) 12,12,11 

11  y*YM0NE*SI(M)-TW0*CL  (M)  >*FM*P(K)/<  <FN*1,0>*CLAT) 

G0  T 0 13 

12  v*Y*<QNE*SL(M)-TWC*CL(M)  )*9(K)«SLAT 

13  L*L*2 

14  M * M ♦ 1 

15  CONTINUE 

CONVERSION  T0  COORDINATE  SYSTEM  SPECIFIED  BY  I TYPE 


3 N E * X 

XsX«CD*Z*SD 

Z=Z«CD-3NE«SD 

T*S'QRT<X*X*Y.Y4Z*2) 

RETURN 

LIBRARY  FUNCTIONS  USED  BY  THIS  SUBROUTINE  aRE  S I n’,  C0S‘, SORT 
SAMPLE  RESULTS 

I OR F( 1965, G, 1,0. 0,0.0,30,0 1 2541. 240,, 563497. 56 40 7,  ) 

I TRF11965, 0,1. 0.0, 180. 0,295. 0,12923, ,9661, *56900, ,59144,) 

I NRF(  1969, 3, 1,0.0, 57. 0,lV5*,  Ci  25092.  ,57297,30145,  ,39638,  ) 

I(, f.R  ( 1965 ,0,2, 63  71 .2,45.0,  >0.  0,21848,  ;-1052,  , 40417  ,^45956, ) 
I ti'  F(  1967, 5, 2, 9,100. 0,31. 0,359. 0,5750., -1226,  .17853.  .18796, ) 


END 


JS164600 

JS104900 

JS185000 

JS185100 

JS185200 

JS165300 

JS185400 

JS185500 

JS185600 

JS185700 

JS185800 

JS185900 

JS186000 

JS186100 

JS136200 

J5186300 

JS186400 

JS186500 

JS186600 

JS166700 

JS166800 

JS186900 

JS187000 

JS187100 

JS187200 

JS187300 

JS1B7400 

JS187500 

JS187600 

JS187700 

JS187800 

JS187900 

JS168000 

JS188100 

JS168200 

JS188300 

JS188400 

JS168500 

JS188600 
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ORCC-RAM  LENGTH 
CNT^Y  POINTS  I GRF 
EXTERNAL  SYMBOLS 

02007110 
080DICT, 
SQRTF 
S INF 
C0SF 

00163  SYMBOLS 

LOAD 

RUN.60.R500 


10ENT 

01266 

00«l3 


IGRf 


j 


47 


